home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 1998 November / IRIX 6.5.2 Base Documentation November 1998.img / usr / share / catman / p_man / cat3 / complib / dlansp.z / dlansp
Text File  |  1998-10-30  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAANNNNSSSSPPPP((((3333FFFF))))                                                          DDDDLLLLAAAANNNNSSSSPPPP((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLANSP - return the value of the one norm, or the Frobenius norm, or the
  10.      infinity norm, or the element of largest absolute value of a real
  11.      symmetric matrix A, supplied in packed form
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      DOUBLE PRECISION FUNCTION DLANSP( NORM, UPLO, N, AP, WORK )
  15.  
  16.          CHARACTER    NORM, UPLO
  17.  
  18.          INTEGER      N
  19.  
  20.          DOUBLE       PRECISION AP( * ), WORK( * )
  21.  
  22. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  23.      DLANSP  returns the value of the one norm,  or the Frobenius norm, or the
  24.      infinity norm,  or the  element of  largest absolute value  of a real
  25.      symmetric matrix A,  supplied in packed form.
  26.  
  27.  
  28. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  29.      DLANSP returns the value
  30.  
  31.         DLANSP = ( max(abs(A(i,j))), NORM = 'M' or 'm'
  32.                  (
  33.                  ( norm1(A),         NORM = '1', 'O' or 'o'
  34.                  (
  35.                  ( normI(A),         NORM = 'I' or 'i'
  36.                  (
  37.                  ( normF(A),         NORM = 'F', 'f', 'E' or 'e'
  38.  
  39.      where  norm1  denotes the  one norm of a matrix (maximum column sum),
  40.      normI  denotes the  infinity norm  of a matrix  (maximum row sum) and
  41.      normF  denotes the  Frobenius norm of a matrix (square root of sum of
  42.      squares).  Note that  max(abs(A(i,j)))  is not a  matrix norm.
  43.  
  44.  
  45. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  46.      NORM    (input) CHARACTER*1
  47.              Specifies the value to be returned in DLANSP as described above.
  48.  
  49.      UPLO    (input) CHARACTER*1
  50.              Specifies whether the upper or lower triangular part of the
  51.              symmetric matrix A is supplied.  = 'U':  Upper triangular part of
  52.              A is supplied
  53.              = 'L':  Lower triangular part of A is supplied
  54.  
  55.      N       (input) INTEGER
  56.              The order of the matrix A.  N >= 0.  When N = 0, DLANSP is set to
  57.              zero.
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDLLLLAAAANNNNSSSSPPPP((((3333FFFF))))                                                          DDDDLLLLAAAANNNNSSSSPPPP((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      AP      (input) DOUBLE PRECISION array, dimension (N*(N+1)/2)
  75.              The upper or lower triangle of the symmetric matrix A, packed
  76.              columnwise in a linear array.  The j-th column of A is stored in
  77.              the array AP as follows:  if UPLO = 'U', AP(i + (j-1)*j/2) =
  78.              A(i,j) for 1<=i<=j; if UPLO = 'L', AP(i + (j-1)*(2n-j)/2) =
  79.              A(i,j) for j<=i<=n.
  80.  
  81.      WORK    (workspace) DOUBLE PRECISION array, dimension (LWORK),
  82.              where LWORK >= N when NORM = 'I' or '1' or 'O'; otherwise, WORK
  83.              is not referenced.
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.